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Nonlinear reaction-diffusion systems are known to exhibit very many novel spatiotemporal 
patterns. Fisher equation is a prototype of diffusive equations. In this contribution we 
investigate the integrability properties of the generalized Fisher type equation to obtain 
physically interesting solutions using Lie symmetry analysis. In particular, we report several 
travelling wave patterns, static patterns and localized structures depending upon the choice 
of the parameters involved. 

■ 

1 Introduction 

■ Nonlinear partial differential equations are frequently used to model a wide variety of phenomena 

ff^ . in physics, chemistry, biology and other fields In such models, when large aggregates 

\ of microstructures consisting of particles, atoms, molecules, defects, dislocations, etc. are able 

' to move and/or interact, the evolution of the concentration of the species can be shown to obey 

Tij- ■ nonlinear diffusion equations of reactive type. These equations play an important role in dis- 

! sipative dynamical systems. Many interesting physical phenomena, such as wall propagation 

in liquid crystals, nerve impulse propagation in nerve fibres, pattern formation in dissipative 
systems, nucleation kinetics and neutron action in the reactor, are closely connected with the 
study of nonlinear diffusion equations. The underlying systems give rise to very many sim- 
ple/complex patterns which are essentially distinct structures on a suitable space-time scale 
^ ' and they arise as collective and cooperative phenomena due to the underlying large number of 

. constituent subsystems. These structures tell us a lot about the dynamics as well as about the 

microscopic behaviour of the underlying systems to some extent. As the interactions among the 
constituents are nonlinear, novel structures which can mimic naturally occurring patterns arise. 
These structures can be stationary or changing with time. 

Generally, in the study of dissipative systems, one of the challenging problems is the selection 
mechanism. That is, one would like to know the kinds of evolving velocity and emerging patterns 
that would be selected in a kinetic process when the system is suddenly quenched into an unstable 
state. Aronson and Weinberger's work on the Fisher type nonlinear diffusion equation j4| has 
shown the existence of distinct selection mechanism, that is the solution u{x, t) of the Fisher 
equation in (1 -|- 1) dimensions, 

Ut = Uxx + U{1 - U), (1) 

converges to a local travelling wave with a definite speed from a wide class of initial data. Further 
it is known that equation has a travelling wave solution called a cline |2| which is nothing 
but a wave travelling in the a;-direction with c > Cmin = 2. However, the first explicit analytic 
form for a cline solution was obtained by Ablowitz and Zeppetella 5 , who showed that an exact 
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propagating wavefront solution (see Fig. ^ is of the form 



u{x, t) = 1 



IH — 1= exp 
\/6 



76^ 



V6 



(2) 



where k is an arbitrary constant. Here the authors made use of the Painleve singularity structure 
analysis of equation ^ to find the exact solution in the year 1979. 
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Figure 1. An exact wavefront solution ^ of the Fisher equation = x — -^j 

There is continuing interest in recent literature j)\ to investigate more general forms of the 
Fisher equation. For instance, there is an interesting generalization of the Fisher equation in 
the description of bacterial colony growth, chemical kinetics and many other natural phenomena 
and it is of the general form |lUj 



du{f, t) 
dt 



(3) 



where D is the diffusion coefficient, A is the nonlocal growth rate, A is the local growth rate, 
G{u, f, t) is the local growth function, V and A are gradient and Laplacian operators respectively. 
As a special case of equation ©, we obtain the generalized Fisher type equation 



ut — An 



m 



1 



■{Vuf - u{l - u) = 0, 



(4) 



where the subscript denotes partial differentiation with respect to time. In the study of popu- 
lation dynamics, u{f,t) refers to the population density at point f at time t. In equation (jlj), 
the linear term modelling the birth rate gives rise to an exponential growth in time while the 
quadratic term that models competition between individuals for food, etc. leads to a stable, 
homogeneous value u = 1 at long times and the diffusion term models the spatial variation of 
the population. This introduces the possibility of spatial pattern formation between the homo- 
geneous regions with u = 1 and u = for appropriate initial conditions. Further, the classical 
Fisher equation (m = 0) occurs in models of population growth [2], neurophysiology [Jj, Brow- 
nian motion [HJ and nuclear reactors 9_. Besides allowing for exact solutions, the m = 2 case 
finds its application in real systems such as the bacterial colony growth ^01 where the square- 
gradient term corresponds to the nonlocal growth occuring at concentration gradients which 
is similar to the nonlinear terms in the Kuramoto-Sivashinsky equation for propagating flame 
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and in the theory of growing interfaces. Moreover, models which admit exact solutions are of 
considerable importance for understanding general behaviour of nonlinear dissipative systems. 
In one dimensional space, such models have received considerable attention. But many realistic 
models are two or three dimensional in nature and in this direction, Brazhnik and Tyson 
considered equation @ in two spatial dimensions and explored five kinds of travelling wave pat- 
terns namely plane, V and y-waves, a separatrix and space oscillating propagating structures. 
All these structures were found when the medium is unbounded and spatially homogeneous. 
Further they show that when the medium is bounded and no flux is allowed through the bound- 
aries, only plane and oscillating waves survive because the frontline of the wave must approach 
the boundary orthogonally. 

In general, obtaining solutions for reaction-diffusion systems is more complex than that for 
pure dispersive systems. For the latter there are several analytical methods like the inverse 
scattering transform method the Hirota method ^21; Backlund transformation method |13j . 
Lie-Backlund symmetries method and so on. On the other hand, for nonlinear diffusive systems, 
no such formal techniques are available to solve them analytically. Very often perturbation 
analysis or numerical techniques are used to treat them. There is therefore an urgent need 
to isolate and identify integrable nonlinear reaction-diffusion systems which can act as model 
systems to deal with more complicated cases. In this connection, symmetry analysis can play a 
very crucial role. 

Consider for example, the well known case of Burgers equation 

Ut = UUxx + UUx, 

where u is the diffusion coefficient. It can be considered to be integrable in the sense that 
it is linearizable: Under the Cole-Hopf transformation u = —uvx/v, it reduces to the linear 
heat equation. It possesses interesting Lie point symmetry structures and infinite number of 
Lie-Backlund symmetries. So, it will be quite interesting to know about other such integrable 
reaction-diffusion equations and the role of symmetries that allows the system to exhibit different 
spatiotemporal patterns and structures which usually possess some kind of symmetry. In this 
direction, the method of Lie groups is the most powerful method to analyse nonlinear partial 
differential equations (PDEs) and hence we make use of it and the singularity structure analysis 
to investigate the integrability properties and hence the dynamics/patterns of We report 
in this paper that the m = 2 case of equation (|1J) possesses infinite dimensional Lie symmetry 
structure, which allows one to linearize it both in (1 + 1) and (2 + 1) dimensions and to obtain 
a large class of exact solutions. We also obtain several exact solutions for the m ^ 2 case. 

The plan of the paper is as follows. In Section 2, we briefly recall some of the important 
reaction-diffusion equations exhibiting novel/complex patterns. Then in Section 3, by carrying 
out the singularity structure analysis, we point out that the PDE Q is free from movable 
critical singular manifolds for the specific value m = 2. More interestingly, we point out that 
the Backlund transformation deduced from the Laurent expansion gives rise to the linearizing 
transformation for this case in a natural way. In Sections 4 and 5, we discuss different underlying 
patterns via symmetry analysis and similarity reductions for the generalized Fisher type equation 
in 1- and 2-spatial dimensions, respectively. Finally we summarize our results in Section 6. 

2 Reaction-diffusion systems and various patterns 

The general form of the nonlinear reaction-diffusion equation is given by 



-^ = V-iDVC) + F{C^,f,t) 
D = dmg{Di,D2,...,Dn), 
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Here C_ represents the population or concentration densities of the species and D and F are, 
in general, nonlinear functions of C_ representing the diffusivity and the reaction kinetics re- 
spectively. In such a case, the dynamics is dominated by the onset of patterns. Inspite of the 
absence of rigorous analytical tools as in the case of soliton systems, combined local analysis 
and numerical investigations on such systems have been found to exhibit a number of important 
spatiotemporal patterns. 

Some of the dominant patterns exhibited by these systems are homogeneous or uniform steady 
states, travelling waves, spiral waves, Turing patterns (rolls, stripes, hexagons, rhombs, etc.), 
localized structures, spatiotemporal chaos and so on. A few of the well known models include 
the following: 

2.1 The Oregonator model 

This model explains the various features of the Belousov-Zhabotinsky reaction and was intro- 
duced by Fields, Koros and Noyes of University of Oregon, USA in 1972. In its simplest version 
it reads as ^3] 



bu2{ui - a) 

- ui) - 



(ui + a) 

U2t = D2^'^U2 + ui - U2- (5) 

Here ui is the concentration of the autocatalytic species HBr02, U2 is the concentration of the 
transition ion catalyst in the oxidised state Ce^"'' and Fe^"*" and rj, a and b are parameters. 
This model is the most popular among the pattern forming chemical reactions. In particular, 
© exhibits 'propagating pulse solutions' that can travel through the system without attenuation. 
Besides, it admits periodic wave trains, target patterns and in two dimensions they generate 
spiral waves. 



2.2 Gierer— Meinhardt model 

It describes possible interaction between an activator a and a rapidly diffusing inhibitor h and 
is of the form jl5j 

(1 + Kaa^) 
ht = DhV'^h + phO^ - + ah, 

where Da and Df, are the two diffusion coeffcients, pa and pb are the removal rates and Ca and dh 
are the basic production terms of the activator and inhibitor respectively. Further Ka corre- 
sponds to the saturation constant. This model is mainly used in the study of the development 
of an organism in biological pattern formation. They are also used to model cell differentiation, 
cell movement, shape changes of cells and tissues and so on. 



2.3 Brusselator model 

Among the various reaction-diffusion type model systems, this is one of the best studied models 
for the formation of chemical patterns theoretically ■ It is based on the chemical reactions 

A — > X, B + X — > y, 2X + Y — > 3X, X — > E, 

where the concentration of the species A, B and E are maintained constant. Thus they form 
the real constant parameters of the system. The evolution of the active species X and Y can be 
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described by 



Xt = A - {B + 1)X + X'^Y + DxV^X, 
Yt = BX- X'^Y + L>y V^y, 



(6) 



after proper rescaling. Here Dx and Dy are diffusion coefficients. This model exhibits hetero- 
geneous patterns through Turing instabiUty. 

2.4 Lotka— Volterra predator-prey model 

Taking into consideration the interaction of two species in which the population of the prey is 
dependent on the predator and vice-versa, the model equations (T7| become 

•S"!* = D\S\xx + ai5i — h\S\S2, 
S2t = D2S2XX — (1282 + b2SiS2- 

Here 5i and ^2 are the population densities of prey and predator. Di and D2 are the diffusivities 
of the two populations, respectively. The parameters oi, 02 are the linear ratio of birth and 
death rates of the individual species while 61 , 62 are the nonlinear decay and growth factors due 
to interaction. 

2.5 FitzHugh— Nagumo nerve conduction model 

The Hodgkin-Huxley model describes the propagation of the electrical impulses along the axonal 
membrane of a nerve fibre. FitzHugh-Nagumo nerve conduction equation JH] is the simplest 
version of the above model and is represented by the following set of equations: 



Here the membrane potential is V{x,t), R corresponds to the lumped refractory variable and 
I{x, t) is the external injected current. The parameters a and b are positive constants while c 
stands for the temperature factor. The above model has been widely used to study various phe- 
nomena in neurophysiology and cardiophysiology. This system exhibits travelling wave pulses [H] . 
In particular, the two-dimensional version of ((7| admits ring wave patterns as well as spiral wave 
patterns for a variety of special initial conditions. 

As mentioned in the introduction, symmetries can play a very important role in determining 
the underlying dynamics of nonlinear systems. Particularly they can help to identify integrable 
cases of the above type of reactive-diffusive systems, if they exist. As an important case study, 
we now investigate integrability and symmetry properties of the generalized Fisher type equa- 
tion (jH). 

3 Singularity structure analysis 

This analysis separates out the m = 2 case for both the (1 -|- 1) and (2 -|- 1) dimensions as the 
only system for which the Fisher equation (jlj is free from movable critical singular manifolds 
satisfying the Painleve property 13._. By locally expanding the solution in the neighbourhood 
of the non-characteristic singular manifold (j){x,t) = 0, (j)x-,4't 7^ in the form of the Laurent 
series JHl 




(7) 
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J+p 
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the possible values of the power of the leading order term are found to be 
» P = -2, 

1 — m 
(in) p = 0. 

For all these leading orders, only for the value m = 2 the solution is free from movable critical 
singular manifolds since for p = —1 the leading order coefficient uq becomes arbitrary besides 
the arbitrary singular manifold (p. In all other cases only one arbitrary function exists for m = 2 
thereby leading to special solutions. 

More interestingly, from the Laurent series expansion if we cut off the series at "constant" 
level term, that is j = —p for the leading order p = 1/(1 — m) = —1, m = 2, one can deduce 
the Backlund transformation that gives rise to the linearizing transformation in a natural way. 
Thus, defining the relation 

u = — + ui, (8) 

we demand that if ui is a solution of equation Q) for the case m = 2, then u is also a solution, 
from which the Backlund transformation is deduced. Now starting from the trivial solution, 
til = of we find that the equations for uq and (p in equation (jH)) are consistent for the 
choice Uq = (j), giving rise to the new solution u = 1. This is nothing but an exact solution of 
equation @. Then with ui = 1 as the new seed solution, one can check from equations satisfied 
by Uq and (j) that 

Uo = -l, (j)t - (pxx - (p + = 0. (9) 
Choosing </> = 1 + X) equation @ can be rewritten as the linear heat equation, 

Xt-Xxx-X = 0. (10) 
Thus the transformation 

^ = 1-T^, (11) 
1 + X 

where x satisfies the linear heat equation UlUj). is the linearizing transformation for equation 
in (1 + 1) dimensions for the choice m = 2 in an automatic way. We note that this is exactly 
the transformation given in ref. :20 in an adhoc way. Here we have given an interpretation 
for the transformation in terms of the Backlund transformation. The same transformation Hll|) 
linearizes equation (3) in (2 + 1) dimensions (for m = 2) as well, where x satisfies the two 
dimensional linear heat equation xt — Xxx — Xyy ~ X = 0- Further equation (jllj) transforms 
equation Q in (3 + 1) dimensions to the 3-dimensional heat equation as well; however, we do 
not study the case further here. 

4 Symmetries and integrability properties 

of (1 + 1) dimensional generalized Fisher equation 

The generalized Fisher equation @ in its (1 + 1) dimensional form reads as 

ut-Uxx-z u + u =0. (12) 

I — u 
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An invariance analysis of equation (|12|) under the infinitesimal transformations 

t + eT{t, X, u), 



X 

u 



X = X + eS,{t, x,u), t — >T 
U = u + e(l){t,x,u), e ^ 1, 



separates out the m = 2 case in that it possesses a nontrivial infinite-dimensional Lie algebra of 
symmetries 



T = a, 



(j) = c{t,x){l - uf 



Here a, b are arbitrary constants and c{t, x) is any solution of the linear heat equation ct—Cxx—c 
0. For all other values of m in equation H12|) one gets trivial translation symmetries 



T = a, 



0. 



In order to obtain solutions of physical importance and corresponding patterns, we make use 
of the method of similarity reductions. This leads to the similarity reduced variables for the 
m = 2 case as 



z = ax — bt, 



1 



+ v{z) + fc{t,x)dt' 



(13) 



Using (|T3|) . equation (|12j) can be reduced for the m = 2 case to the similarity reduced ordinary 
differential equation (ODE) 

a^v" + bv' + v = Q 

whose general solution is 



"1-1,2 



-b ± - 4a2 
2^2 ^ 



where Ii and I2 are integration constants thereby leading to 



b^ - 4a^ > 0; 



62 _ 4^2 ^ Q 



a + I^e^riax-bt) _,_ J^^m2(ax-bt) _^ J c{t,x)dt' 
a 

a + eP(''^-^*) (/i + /2(ax - bt)) + / c(t, x)(it ' 

a 

a + eP(ax-W) [I-^^ cos q{ax - bt) + I2 sin q{ax - bt)) + J c{t,x)dt' 



9 - < 



with p = —b/2a?, q = V ^a? — 62/2a^, as the solution to the original PDE (|12j) . Here the 
similarity reduced variable (|13() is nothing but the linearizing transformation (jllj) . 

Proceeding in a similar fashion for all the other (nonintegrable) cases (m 7^ 2), the similarity 
variables z = ax — bt and u = w(z) reduce equation 1)12(1 to the ODE 



aw— ma v + bvv — (1 — v)v 



0, 



1 — w. 



(14) 



which is in general nonintegrable except for m = and b/a = 5/-v/6- This special choice leads 
to the cline solution obtained by Ablowitz and Zeppetella 5 . In the static case (6 = 0), one 
obtains elliptic function solutions. Besides, a particular solitary wave solution 



(3 - 2m) 
(2 - 2m) 



sechM/2-|J^ 



m < 1 
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with I2 as the second integration constant, which is a hmiting case of a elhptic function solution, 
is also obtained (refer Fig. \^ . In the general equation (|14|) is of nonintegrable nature, 

we make use of numerical techniques to study the underlying dynamics. Here we obtain typical 
periodic wave trains for 6/a = which is in accordance with the fact that reaction-diffusion 
systems exhibiting limit cycle motion in the absence of diffusion exhibits travelling wave patterns 
(Fig. EKjb). For h/a = 1, we get a propagating pulse (Fig. |3t) and the corresponding phase 
portrait {v — v') shows a stable spiral equilibrium point (Fig. Eli)- On increasing the value of 
6/a, that is, at 6/a > 2 (6 = 2.041) [5:, the system supports a travelling wave front (Fig. and 
the trajectories in the phase plane {v — v') correspond to a stable node (Fig. I^f). 




-20 20 

X 

Figure 2. A static solitary wave pulse for to = 1/2 of the generalized Fisher equation (|12|l . 



5 The (2 + 1) dimensional generalized Fisher equation 

Extending a similar analysis to the (2 + 1) dimensional case of the generalized Fisher equation 

7T}j / 9 9 \ 9 / \ 

Ut - Uxx - Uyy - Y—^\Ux + Uy) - U + U =0, (15) 

one finds that the invariance analysis of equation H15|) under the infinitesimal transformation 
singles out the special value m = 2 for which the Lie point symmetries are 

T = a, ^ = 63^ + 64, r/ = -63X + d4, (j) = c{t,x,y){l - uf , 

where r] is the infinitesimal symmetry associated with the variable y, c(t, x, y) is the solution of 
the two dimensional linear heat equation q — Cxx — Cyy — c = and 63, 64 and are arbitrary 
constants. But for all other choices of m (7^ 2) we get 

T = a, ^ = 631/ + 64, r/ = -63X + ^4, (t) = 0. 

In a similar fashion as that for the (1 + 1) dimensional case, the similatiry variables for the 
m = 2 case 

^3 / 2 , 2^ , , , . a • -1 / di-bsx \ 

zi = —{x +y) + biy-diX, 2:2 = -t - — sm ' ' 



w{zi,Z2) + f c{t,x,y)dt 



(16) 
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reduce the PDE ^ to 



w 



22 



+ 2b3W,, + (263Z1 + bl + dfj w,,,, + 



aw 



2h3Zi + bl + dl 



+ w — a = Q. 



Here too one can obtain the hnear heat equation 

Xt - Xxx - Xto - X = 0, 



w{zi, Z2) + / c{t,x,y)dt 



X 



from the similarity form ()16|). Such a transformation can be interpreted as the Unearizing 
transformation from a group theoretical point of view. 

Carrying out a Lie symmetry analysis for equation (\n\i also, one can obtain the new similarity 
variables 



( '^1^2 

C = zi, w = a + ^ 



1 



/(C) + - / d2{zi,Z2)ey =3-;dz2 

C3 



zi = 2bszi + b\ + dl, Z2 = Z2, &3, d.4 / 0, 
where / satisfies the linear second order ODE of the form 
ef" + Cf' + {A + BC)f = 0, 



A = (aci/263C3) 



B = {l + ci/c3)/4b 



3i 



(18) 
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with prime denoting differentiation w.r.t. (. Thus the solution to the original PDE reads as 
u = 1 — a 



=1 



^2 



a + eV-3~-/ I I^Zi [2^/BFlj + I2Z2 [2y/mi 

1 



C3 



eV'^a 'dz2 + / c{t,x,y)dt 



(19) 



In the limit 64 = (i4 = ci = the system is found to exhibit circularly symmetric structures 
given in Fig. ^ 



u 




-10 



10 10 



Figure 4. Circularly symmetric patterns of the (2 + 1) dimensional generalized Fisher equation IjlSfl for 
m = 2. 

More interestingly, in the special case 63 = 0, ^4 = the system exhibits propagating wave 
structures and the corresponding forms are 



u = < 



1 — a< a + exp 



k I —X — t 
M 



Ii cos (Vfci 05641/) + /2sin ( v/ci 05642/ 



C5 



+ J c{t,x,y)dt\ , ki<0, 



1 — a< a + exp 



k I —X — t 
64 



C5 



+ J c{t,x,y)dt 



ki > 0, 



1 — a< a + exp 



k I —X — t 



+ 



I 



C3{Z1,Z2) ^kz2 



e^^^dzi 



C5 



+ j c{t,x,y)dt 



h =0, 



where the parameter ki = tt^ 



1 



(20) 



with k = — C2/C5 and zi = b^y, Z2 = f^x — t 



Here C2, C4, C5 are arbitrary constants of integration. Equation (|20|) . in particular exhibits 
the five classes of bounded travelling wave solutions reported by Brazhnik and Tyson ^ for 
certain choice of the parameters involved along with the specific assumptions of the functions 
C3{zi,Z2) = and c{t,x,y) = 0. The corresponding solutions are given below. 
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Among the classes of solutions, the simplest travelling wave solution (Fig. EK) 

1 



1 



1 + A exp 



-k 



b4' 



t]± 



can be constructed by assuming either Ii 
pattern (Fig. [SJa) 

1 



= or /2 



1 



1 + A exp 



-k 



b4' 



t 



cosh [y/kicc^biy) 



fci > 



0. For Ii = I2 (/ 0), we obtain a V-wave 



ki > 0. 



Again the case I2 
represented by 

u = l 



and A;i < leads to a wave front oscillating in space (Fig. [St) and is 



1 + A exp 



t 



I cos {Vklc^h4_y)\ 



But when /i / and /2 = we get a separatrix (Fig. Eli) 

1 



1 



1 + A\y\ exp 
Finally for positive ki and Ji 
n = 1 



-/2 the Y-wave solution (Fig. becomes 
1 



1 + A exp 



t 



sinh (Vhc^b^y) 



In each of the above solutions ^ is a positive constant. It is a well known fact about Fisher 
equation is that it forms a basis for many nonlinear models of different nature. As a result, 
the above solutions are reminiscent of patterns from different fields. In particular, V-waves are 
characterized in the framework of geometrical crystal growth related models |2j and in excitable 
media |22j while space oscillating fronts are relevant to cellular flame structures and patterns in 
chemical reaction diffusion systems |23j . Further it has been shown in |21] with a geometrical 
model that excitable media can support space-oscillating fronts. Several static structures can 
also be obtained as limiting cases of the above solutions (|T9|l and (|20j) . 

Finally, a similar analysis for the nonintegrable (m 7^ 2) case yields static patterns/structures 
in (x, y) variables. Here one has to look for certain special solutions due to its nonintegrable 
nature. That is, for 63 = and ^4 = with the similarity variables zi = 647/, Z2 = f^x — t, 
u = w{zi, Z2), the reduced ODE reads as 



ci/' + /(l-/) = 0, 



^ = ( 12^1 + bW 



(Y = d/dC, 



with C = — ci (^-^^x — + C2biy and w = /(C)) giving rise to plane wave structures. For 63 = 0, 

the similarity variables zi = d^x — b^y, Z2 = ax — b^t and u = w{zi, Z2) reduces the PDF to an 
ODE 



Ji 



+ = 0, 



0' = d/dC 



with fi = l — f,A = a? {c\d\ + C264), B = —^464(01 + C2) and Q = ac2{diX — b^y) — di{ci + 
C2)(ax — 64*), w = f{C)- Then the system is found to possess elliptic function solutions including 
the limiting case of the solitary pulse for certain choices of the constants involved. 
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id) (e) 

Figure 5. Five interesting classes of propagating wave patterns as obtained in ref. 'E', which follow from 
equations H2(J|) : (a) travelling waves; (b) V-waves; (c) oscillating front; (d) separatrix solution; (e) Y-waves 
withe = -fc(^x-i). 

6 Conclusion 

Our studies on the integrability /symmetry properties of the the generalized Fisher type nonlinear 
reaction-diffusion equation show that the system under consideration possesses interesting Lie 
point symmetries that could form infinite dimensional Lie algebra for the particular choice of the 
system parameter m = 2, thereby exhibiting various interesting patterns and dynamics. Besides, 
the singularity structure analysis singles out the m = 2 case as the only system parameter for 
which the generalized Fisher type equation is free from movable critical singular manifolds. The 
generalized Fisher equation is found to possess a large number of interesting wave patterns. It 
will be of interest to consider other physically interesting reaction-diffusion systems from the 
Lie symmetry point of view and to study the underlying patterns. 
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